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Abstract. We exhibit an infinite family of triplets of mutually 
unbiased bases (MUBs) in dimension 6. These triplets involve the 
Fourier family of Hadamard matrices, F(a,b). However, in the 
main result of the paper we also prove that for any values of the 
parameters (a, 5), the standard basis and F{a,b) cannot be extended 
to a MUB-quartet. The main novelty lies in the method of proof 
which may successfully be applied in the future to prove that the 
maximal number of MUBs in dimension 6 is three. 
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1. Introduction 

The notion of mutuaUy unbiased bases (MUBs) emerged in the hter- 
ature of quantum mechanics in 1960 in the works of Schwinger [2H] • It 
now constitutes a basic concept of Quantum Information Theory and 
plays an essential role in quantum-tomography [20l |32] , quantum crip- 
tography ^ElETj, the mean king problem pLj as well as in constructions 
of teleportation and dense coding schemes [3T] . 

Recall that two orthonormal bases of C^, A = {ei,...,ed} and 
B = {fi,...,frf} are said to be unbiased if, for every I < j,k < d, 

|(e,-,ffc)| = —=. A set Bo,...Bm of orthonormal bases is said to be 
Vd 

mutually unbiased if any two of them are unbiased. It is well-known 
(see e.g. [2], El [lEl [19], [231 132]) that the number of mutually unbiased 
bases in cannot exceed d+1. It is also known that d + 1 such bases 
can be constructed if the dimension d is a prime or a prime power (see 
e.g. [21 [121 [131 [H [201 [211 [32]). Apart from this, very little is known 
except for the fact that there are always p+1 mutually unbiased bases 
in where p is the smallest prime divisor of d. Thus, the first case 
where the largest number of mutually unbiased bases is unknown is 
d = 6: 
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Problem 1.1. 

What is the maximal number of pairwise mutually unbiased bases in 
C^? 

Although this famous open problem has received considerable atten- 
tion over the past few years ([5], [9l HOl |26l |29|), it remains wide open. 
Since 6 = 2 x 3, we know that there are at least 3 mutually unbiased 
bases in C^, but so far tentative numerical evidence [HI [IDl IHl E3] sug- 
gests that there are no more than 3, a fact apparently first conjectured 
by Zauner [53]. 

One reason for the slow progress is that mutually unbiased bases are 
naturally related to complex Hadamard matrices (and the classification 
of such matrices in dimension 6 seems to be very difficult). Indeed, 
if the bases Bq, . . . ,Bm are mutually unbiased we may identify each 

Bi = {e^i \ . . . , ef^} with the unitary matrix Uj = 



i. e. the k-th column consists of the coordinates of the k-th vector of Bj 
in the bases Bq. (Throughout the paper the scalar product (., .) of C'^ 
is linear in the first variable and conjugate-linear in the second). With 
this convention, Ui = Id the identity matrix and all other matrices are 
unitary and have entries of modulus l/\fd. Such matrices are called 
complex Hadamard matrices. It is clear that the existence of a family of 
mutually unbiased bases Bq, . . . , Bm is thus equivalent to the existence 
of a family of complex Hadamard matrices Hi, ... , Hm such that for all 
^ ^ j k < m, H*Hk is again a complex Hadamard matrix. In such 
a case we will say that these complex Hadamard matrices are mutually 
unbiased. 

A complete classification of complex Hadamard matrices is only 
available up to dimension 5 (see [IH])- The classification in dimen- 
sion 6 is still out of reach despite recent efforts [3], [211 EH] • This is one 
of the reasons for Problem II. II to be difficult. 

A natural question that arises in this context is that given two un- 
biased orthonormal bases, does there always exist a third orthonormal 
basis that is unbiased to the first two? Or, equivalently, given a com- 
plex Hadamard matrix H, does there always exist another one G that 
is unbiased to H7 The answer is negative in such generality. It was 
recently proved in [TUj that for the matrix Sq (cf. [3DJ for the nota- 
tion) there exists no complex Hadamard matrix unbiased to it. A less 
restrictive question is the following: 

Problem 1.2. 

Given a Hadamard matrix H , does there always exist some unbiased 
vectors to H? 

At this stage, it may be worth recalling some invariants of this prob- 
lem. Assume Hi, . . . , Hm are mutually unbiased n x n Hadamard 
matrices, and let D,Di, . . . , -Dm be unitary n x n diagonal matrices 
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P, Pi, . . . , Pjn henxn permutation matrices. Then DPHiPiDi are still 
Hadamard matrices and are still mutually unbiased. We will say that 
DPH\*^ PiDi is equivalent to where the superscript (*) is a choice 
(the same for all matrices) between complex conjugation or nothing. 
Note that the effect of this operation is that we may assume that the 
first row of each Hadamard matrix is . . . , 1] and that the first 

column of Hi is . . . , 1]. We may further order the remaining 

rows and columns of Hi. 

Further, note that if an n-dimensional vector 

^■^^ y _ ^—(1 g2j7r(/>i g2i7r</>2 ^2m<f>z ^2i-n<f>A ^2in<)>5^ 

is unbiased to the standard basis and to the n — 1 first columns of a 
Hadamard matrix, then it is automatically unbiased to the last one. 
Hence, we have n — 1 unbiased-criteria to be satisfied for the n — 1 
parameters (pj. In generic situations we therefore expect a finite number 
of solutions to arise. We know of a non-generic example (in dimension 
4) where infinitely many unbiased vectors arise, but of no examples 
where the number of such vectors is zero. 

Moving back to the 6-dimensional case we note the significance of 
Problem II. 2[ If for a certain complex Hadamard matrix H the number 
of unbiased vectors is less than 30, then the MUB-pair {Id, H} can 
obviously not be extended to a full set of 7-MUBs (because we would 
need at least 30 vectors to form another 5 bases). 

Further, it is easy to see that a vector u of the form ([T]) is unbiased 
to the columns of Id and of H if and only if the mapping 

6 

(2) 7^(01,02,03,04,05) Ku,hj)|, 

i=i 

(where hj denote the columns of H) has a global maximum at the point 

(01,02,03,04,05) in 

Therefore, a natural way to search numerically for unbiased vectors u 
is to start from a random point of the parameter space (0i, 02, 03, 04, 05) 
in and find local maxima of 7i defined in It is plausible to expect 
that if we run our numerical search many times we will find most, or 
indeed all, unbiased vectors u in this manner (as well as finding possible 
other local maxima which we simply discard). There is no guarantee, 
of course, and we will need to back up our numerical evidence with 
rigorous mathematical statements. 

In most of this paper, we will focus on H belonging to the "Fourier 
family". Let us recall (cf. [30]) that this is the two-parameter family 
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of complex Hadamard matrices -F(a, b) defined by 

1 ■ 

-ux 

-1 ' 

LUX 

where x = e^™°- ^ y = e^'^'"^ and u = e^*''/^. 

Note that -F(0, 0) = jFg is the standard Fourier matrix. We may 
thus see the task of finding unbiased vectors to the standard and the 
F{a, b) bases as a perturbation of the so-called Pauli Problem. Recall 
that Pauli asked whether a function / in L'^{M.'^) is uniquely determined 
by its modulus |/| and the modulus of its Fourier Transform |/| (see 
e.g. [21] or [15] for some results and further references). The discrete 
analogue of this problem, i.e. the problem of finding finite sequences 
of complex numbers of modulus one (aj) such that their Fourier trans- 
forms have also modulus one (such sequences are called hiunimodular 
sequences) has been considered e.g. in [HIH]. Our problem can thus be 
seen as a perturbation of the discrete case in dimension 6. 

For the particular case of H = J-'q, a full analytical solution of Prob- 
lem 11.21 is actually known [3 [17] : there are exactly 48 vectors, nor- 
malized as in ([T]), that are unbiased with respect to {/c/, jFg} and one 
can form 16 different orthonormal bases Ci, . . . Ciq out of them. How- 
ever, no pair of bases (Cj, Cj) are unbiased with respect to each other, 
which means that no triplet {/, jFg, C} can be extended to a mutually 
unbiased-quartet {J,jFg,C, i^} (see [I71[TU|). What happens if we set 
H = FQ{a,b) for some generic values a, 6? We heuristically expected 
that in such a case significantly less than 48 unbiased vectors u should 
arise. We also expected that only in exceptional cases should there ex- 
ist a basis C built from these unbiased vectors. These heuristics turned 
out to be falseQ 

We ran the numeric search of finding local maxima of expression 
(|2]) for several values of a, b. The results were both surprising and 
overwhelmingly convincing. 

Numerical Evidence 1.3. 

For any values of the parameters {a,b), the number of vectors unbi- 
ased to the identity matrix and F{a,b) is 48. For generic values {a,b) 
there are 8 different orthonormal bases Ci{a, b), . . . C^{a, b) that can be 
formed out of these 48 vectors. For some exceptional values of (a, b) 
there are more such bases: for (a, 6) = (0,0) there are 16, while for 



^The authors are grateful for W. Bruzda for making the first computer search, 
with a method other than maximizing expression ([2]), which indicated that the 
number of unbiased vectors is more than 40 on average. 
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(a, 6) = (1/6,0) there are 70 such base^. However, no mutually unbi- 
ased triplet of the form {Id, F{a,b),C) can be extended by a further 
basis to form a mutually unbiased quartet {Id, F{a, b), C, D). 

We will back up most of these numerical data by rigorous analytic 
results in subsequent sectionsH 

In Section [2] we construct an infinite family of MUB-triplets in ana- 
lytic form involving the Fourier family of Hadamard matrices F{a, b), 
with a = (and some restrictions on b). We have recently been in- 
formed by G. Zauner that his work ^3J also includes an infinite family 
of MUB-triplets, although the formulas are not made explicit. As the 
beautiful construction of is scarcely known and it is originally writ- 
ten in German we decided to provide an English version of it in the 
Appendix of ^22j. We will show, however, that Zauner's family is not 
equivalent to ours. 

One may think that the emergence of an infinite family of MUB- 
triplets is a major step towards finding a MUB-quartet in dimension 
6. On the contrary, we prove the following: 

Theorem 1.4. 

None of the pairs {id, F{a, 6)) of mutually unbiased orthonormal bases 
can be extended to a quartet {^Id,F{a,b),C,D^ of mutually unbiased 
orthonormal bases. 

While this can be disappointing for some, we believe that this is a 
breakthrough result of the paper in that the method we apply here may 
later be generalized to settle Problem \l . 1\ and prove that the maximal 
number of mutually unbiased orthonormal bases in dimension 6 is three. 

The paper is organized as follows: Section [2] is devoted to charac- 
terizing vectors unbiased to the standard basis and F{a,b), and the 
construction of an infinite one-parameter family of MUB triplets. In 
Section [3] we prove Theorem 11.41 Finally, in Section H] we attempt to 
offer some general theoretical reasons behind our results, other than 
just the sheer numbers and formulae. 

2. An infinite family of MUB-triplets involving F{a,b) 

In this section we first describe a reduced system of equations for 
any vector u unbiased to the bases A = Id and B = F{a, b). However, 
we can only obtain some particular solutions in closed analytic form 

^The exact number 70 is given in [10] 

"'The authors of [TO] have used the technique of Grobner bases to prove that 
the number of unbiased vectors is indeed 48 for several (but finitely many) tested 
values of {a,b). In fact, in [10] several members of all known Hadamard families 
are tested, not only the Fourier family. However, the techniques of this paper have 
the advantage that they enable us to reach rigorous conclusions about the whole 
family F{a, b) and not just the tested values of the parameters. 
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in the special case a = 0. Nevertheless, these explicit formulae give 
a rigorous proof of the existence of an infinite 1-parameter family of 
MUB-triplets involving A = Id and B = F{0,b). Recall, that the 
numerical evidence actually suggests existence of MUB-triplets for all 
a,b, with A = Id and B = F{a,b), i.e. a two-parameter family. We 
cannot give a rigorous argument in such generality. 

2.1. A reduced system of equations for unbiased vectors. 

We begin with a useful lemma about vectors unbiased with the 
Fourier basis in C^: 

Lemma 2.1. 

Let u = e^*'^/^ and let a,P,-f E C. Then 

(\a + (3 + 7|2 =6 

(4) < |a + + 7c^T = 6 

[|a + Pu^ + 7u;p = 6 

if and only if 

r|a|2 + |/3|2 + |7|2 = 6 

I a/3 + /37 + 7a = ■ 

Proof. One easily sees that (jl]) is equivalent to 

a|2+ |7|2 + 2Re(a/5 + /37 + 7a) = 6 
a|2 + + |7|2 + 2Re{Lj^ap_+uj^p^ + u^-fa) = 6. 
ap + |/9p + |7p + 2Re {u a f3 + uj(3'y + uj-f a) = 6 

Then, using the fact that l + tu + tu^ = and adding all three equations, 
we see that this is equivalent to 

I«l' + l/3|' + l7l' = 6 
Re(a/? + /37 + 7a) = 0. 
Re[ua(3 + uP'j + wyd) = 

We conclude by noticing that R.e{z) = and R.e{ujz) = if and only if 
^ = 0. □ 

Let us now assume that u G is a unit-norm vector that is unbiased 
to the standard basis: 

u — Ci, C2, C3, C4, C5), |ci| — IC2I — IC3I — IC4I — Icsl — 1, 

where the conjugate signs are introduced for later convenience of cal- 
culations. 

Recalling the notation x = e^*™, y = e'^^'^^, the vector u is further 
unbiased with respect to the generalized Fourier basis F{a,b), if and 
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while applying it to Equations (!6l).(l8l).(fT0l) we obtain the following: 

(13) |l-C3p + |x(-Ci+C4)|'+|y(c2-C5)|2 = 6 
(1 - C3)X(C4 - Ci) + X(C4 - Ci)y{c2 - C5) 

(14) +y(c2-C5)(l-C3) = 0. 

But, using the fact that ci, . . . ,05 are all of modulus 1, we see that 
Equation (ITTl) is equivalent to 

(15) Re(c3 + Cici + csc^) = 0. 

Similarly, as |x| = \y\ = 1, ffT^ reads |1 — C3P + |c4 — Cip + |c2 — C5P = 6 
which also reduces to ffT^ . 

We have thus proved the following lemma: 

Lemma 2.2. 

A vector u = ^(1,01,02,03,04,05)) is unbiased to both the standard 
basis and the generalized Fourier basis F{a, h) if and only if all cj have 



absolute value 1, and the following conditions are 


fulfilled 




Re (03 + O1O4 + O5O2) = 


m 




(1 + 03)(0i + O4) + (Oi + 04)(02 + O5) 




< 


+ (02 + 05)(l + 03) =0 


m 




(1 - 03)a;(o4 - oi) + a;(o4 - Ci)y{c2 - 05) 






+1/(C2 - C5)(l - O3) = 


(14) 



Remark 2.3. The most natural way to prove all statements indicated 
by Numerical Evidence 11.31 would be to solve this system of equations 
and show that there are exactly 48 solution vectors for any choice of 
a, b. This would provide an exhaustive list of MUB-triplets involving 
the generalized Fourier matrices F{a,b). Unfortunately, we have been 
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unable to fulfill this task in this generality and we are not able to give 
all solutions in closed analytic form. 

2.2. An infinite family of MUB-triplets involving F(0, b). In or- 
der to obtain analytic formulae for some of the arising MUB-triplets we 
need to restrict our attention to the case a = 0. Even in this case the 
calculations are rather long and cumbersome, and not very instructive. 
The full details are presented in the Appendix of [22], and we only in- 
clude the final result here. Note that all emerging formulae are explicit 
so that the correctness of the result can be checked (most conveniently 
by computer algebra) without going through the detailed calculations. 

Theorem 2.4. Assume that | arcsin ^ l£ t < ^ — ^arcsin^. Intro- 
duce the following variables: 

, , , V2 + cos 2t cos2t / cos^2t 
(16) ^ = arccos ^ , 03 = ^ h^ll ^ — 



1/2 



(17) 



= arccos 



V3 + sin^ 2^ + 3^9 sin^ 2t - 5 



8 sin 2t 

Then define ip and (p by the equations 

— cos^ ip COS t + cos(/? + ip) sin ip sin t 
sin (3 sin 2t ' 

sin cos cos t — sin(/5 -|- if)) sin if) sin t 



(18) 

(19) 
and 
(20) 

(21) 
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sin f3 sin 2t 



— sin t sin^ ip + cos ip cos t sin(/3 + ip) 

sin (3 sin 2t 
cos ip cos(/3 + ip) cos t — cos ip sin t sin ip 
sin (3 sin 2t 



Write r] = e**, u = e^"^ , ^ = e"^ , and b = -^p. Finally define C{t) to be 
the orthonormal basis given by columns of the matrix 
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r/ien t/ie standard basis, the generalized Fourier basis F(0,b{t)) and 
the basis given by C{t) are mutually unbiased. 
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Remark 2.5. This theorem exhibits an infinite family of MUB triplets 
in terms of a parameter t. Each member of the family contains the 
standard basis and one member of the family F{0, b{t)). However, the 
dependence of t on 6 is only implicit and seems unsolvable in closed 
form. Note also that b{t) does not take the value 0. 

We do not claim that we have found all MUB-triplets containing 
the standard basis and F(0,b{t)), but only one such triplet. Actually, 
Numerical Evidence 11.31 shows that there exist other solutions, but we 
have been unable to describe them all analitically. 

Note also, that the family of MUB-triplets above is different from 
the one presented in [33j. This fact is shown in the Appendix of \22\ . 



Remark 2.6. It is natural to ask whether C{t) provides a new family 
of complex Hadamard matrices of order 6. This is not the case as C{t) 
also belongs to the the generalized Fourier family F{a,b). This can 
easily be seen by dephasing the first column and properly reordering 
the remaining ones. 

3. No QUARTET OF MUTUALLY UNBIASED BASES INVOLVING THE 

IDENTITY AND F{a,b) 



In this section we prove Theorem II. 4[ i.e. the non-existence of quar- 
tets of mutually unbiased bases of the form {Id, F{a, b), C, D) for any 
values of a, b. This will be done via a discretization scheme and an 
exhaustive computer search after establishing proper estimates of the 
error terms. We believe that the method can be generalized in the 
future to prove that the maximal number of MUBs in dimension 6 is 
three. 



Proof of Theorem L4_. Let us briefiy describe the basic idea and turn to 
the details later. The proof proceeds by contradiction: assume there 
exists a MUB-quartet {Id, F{a,b),C, D). First, as described in the 
introduction, we may take advantage of the equivalence relations of 
Hadamard matrices to reduce the range of parameters. A priori, (a, b) 
is any point in the square [0, 1)^. But, due the equivalences described 
in [3] we can assume without loss of generality that (a, b) lies in the 
triangle T with vertices (0,0), (1/6,0) and (1/6, 1/12) which is a fun- 
damental region (see for details). 

Next, e^*'^°, e^*'^'' and all entries of \/6C and \/6D are unimodular 
complex numbers. We will thus approximate them by iV-th roots of 
unity and replace the matrix F{a, b) by a matrix F{d, b) and a/6C, 
by matrices V^C, VGD with entries exclusively A^-th roots of 
unity. Of course, in doing so, we will destroy the main features of C, D: 
namely, C and D are neither unitary nor unbiased to F{d, b) (or to each 
other) anymore. However, if A^ is large enough, then {Id, F{d, b), C, D) 
will still approximately be a MUB-quartet. Moreover, the bounds in 
these approximations can be precisely controlled. It turns out that 
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if is large enough, an exhaustive computer search shows that no 
quartet of matrices satisfies the prescribed bounds. This means that 
the hypothetical quartet {Id, F{a,b),C, D) cannot exist. The code of 
the computer algorithm and the full documentation of the results are 
available at the web-page [M] . The running time of the code was about 
6 hours on a computer with a 3,2 GHz CPU. 

Let us now describe the details. Let be an integer. We parti- 
tion the interval [0, 1) into N sub-intervals • • • , I^^^i of equal 
length. I.e. if^ = [j/N, (j + l)/iV), and denote by rf^ = (j + 1/2) /N 
the midpoint of /j^''. Now, if a (resp. b) fall in some interval jj^'' (resp. 

Im^), we then replace a by a = rj^-* (resp. 6 by 6 = rm^). When doing 
so we must keep in mind that the actual value of a (resp. b) can lie 
anywhere in the interval Ij'^'^ (resp. I^^). 

Next, recall that there is no loss of generality in assuming that 
all vectors in all appearing bases have first coordinate 1 / ^/6. All 
other entries have modulus 1/a/6. This has the consequence that 
all appearing scalar products throughout this section have the form 

(u,v) = |(l + ELie^^^^'^'=-^'=^)- 

Denote the entries of \/6C, ^/6D as c^j = e^™'^^'^ and d^j = e'^™^'''^ 
with 'jkjyPkj G [0,1), and the indexing being set as < k,j < 5. 
Thus, for A; = we have •ykj = Pkj = 0, and for /c > 1 each 'jkjyPkj 
falls into some interval ■* , where A^' is another integer (for clarity 
of notation the dependence of i on 'yk,j,Pk,j has been dropped). We 
define VQC, VQD by replacing these entries by r^^ It turns out that 
we can take A^' smaller than A^, which saves a lot of computing time. 
Actually, our search was carried out with A^ = 180 and A^' = 19. 

Finally, the algorithm runs in two steps. In the first one, we seek all 
vectors u of the form 

(22) = ^(1, e2-(^-^+V2)/^', . . . , ^2^.U,+l/2)/N'^ 

V6 

that are "almost" unbiased to F{d, b). These vectors are the candidates 
for the columns of C,D. The second step then consists of constructing 
"almost" orthonormal bases C, D out of these vectors, and checking 
whether those can possibly be "almost" unbiased to each other. Of 
course, all the "almost" terms above need to be properly quantified. 

Let us now turn to the error term. We want to approximate a column 
vector of C or 

(^23) u = ——(1 e^*'^'^^ g2i7r</>2 g2i7r03 ^2in<f>4 g2i7r05 

a/6 ' 

by a vector u of the form where each 0^ has been approximated 
by some {jk + l/2)/N' = ^ ■ We must keep in mind that 0^ can 
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lie anywhere in the interval [jk/N',{jk + 1)/N'). Let us denote by 
fo, . . . , f5 the columns of F(a, b) and by fo, . . . , those of F(a, b). By 
construction u and fo, . . . , fs have the following property: there exist 
numbers (pk in [jk/N', {jk + l)/N'), and numbers a, 6 in [a — j^,d+ j^), 
\b ~ b + such that the corresponding vectors u (as in fpHj) ) and 
fo,...f5 (as in the columns of ([3])) are unbiased to each other, i.e. 
|(u, ffc)| = I/VG. For a fixed pair of discretization parameters N,N' 
and a fixed pair of a, b the vectors u of the form ( |22l) which have 
the above property will be called quasi-unbiased to F{a,b), and their 
set will be denoted by FUB'^J'^,. Our first aim is to find the vectors 
belonging to FUB'^^^, out of all vectors u of the form ([22]) • Despite 

having N'^ possible vectors u, we will see that the set FUB'^^^, will be 
of reasonably small size. 

We will need the following lemma which we will refer to as the trivial 
error bound. 

Lemma 3.1. Let Ik and Jk ^ k < 5) be closed intervals (possibly 
degenerate) contained in [0, 1]. Let L^ andT^ denote the lengths, while 
uik and Sk the midpoints of the intervals Ik and Jk, respectively. 

Consider the midpoint sum S = ^(1 + ^^=16^*'^''™''="'''=^). The follow- 
ing two statements hold: 

• if it is possible to select points (pk o-nd ipk from the intervals Ik 
and Jk, such that ||(1 + ^^^^ e^^^^'^'^"'^''))! = then 

(24) lY.{Lk + Tk)> 

k=l 

• if it is possible to select points (pk CLnd ipk from the intervals Ik 
and Jk, such that |(1 + Y!'k=i e^^'^^'^*-^'^)) = 0, then 

5 

(25) ^^(L, + T,)>|^|, 

k=l 

Proof. Let us introduce the "error function" E{x,y) = 5 — |(1 + 
e^*'"^^*"^''^), where Xk,yk are in Ik and Jk, respectively. Note 
that \{mk — Sk) — {xk — yk)\ < ^i^k + Tk) for each 1 < k < 5. The 
trivial estimate \e'^i^('^k-sk) _Q'2in(x„-yk)^ < vr(Lfc + Tfc) yields \E{x,y)\ < 

f T.l=i{Lk+Tk). Therefore, the values of the function |(1+ELi e^^'^^'^^"^^)) 
stay within a disk of radius | Yfk=ii'^k + Tk) around S. In the first 
statement of the lemma we thus conclude that the distance of S from 
the circle of radius l/y/Q (centered at the origin) is not greater than 
f Ylk=ii-^k + Tk). In the second statement we conclude that the dis- 
tance of 5* from the origin is not greater than | Ylt=iiTk + Tk). These 
are equivalent to and fl25l) . □ 
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Let us apply this lemma to our particular case. The last 5 coordinates 
of u represent intervals of length 1/A^'. For j = 0,2,4 the columns tj 
do not contain the parameters a, b so that all coordinates are known 
exactly, which means that the corresponding intervals are degenerate 
(of length zero). For j = 1,3,5 the columns fj contain parameters 
at four different coordinates, each representing intervals of the type 
[d — a + and [b — b + 2^), of length 1/N. Therefore Lemma 
13.11 yields 



(26) 



(27) 



< 7=0,2,4. 



Stt An . , „ ^ 
< + 7 = 1,3,5. 



However, these bounds turn out to be too crude, and we will need the 
following improved error bound. The technical lemma below establishes 
the simple fact that "maximal error always occurs at the endpoints of 
the intervals". 

Lemma 3.2. Let Ik and Jk ^ k < 5) be closed intervals (possi- 
bly degenerate) contained in [0,1]. Let Lk and Tk denote the lengths, 
rrik and Sk the midpoints, and i^,ik O'^'^d j^jjk endpoints of the 
intervals Ik and Jk, respectively. Assume that Y^k=ii^k + Tk) < ^. 

Consider the midpoint sum S = |(1 + X]fe=i 6^*'^*-™'="'*'=''), and all the 

32 endpoint sums = |(1 + Y^k=i e^*'^'-*^'' ''^) (where e denotes any 
vector of ± signs; note that —tk appears at the upper index of jk). The 
following two statements hold: 

• if it is possible to select points (pk and ipk from the intervals Ik 
and Jk, such that ||(1 + ELi = 7g' ^^^n 

(28) max{|5-5e|} > -^1, 

v6 

• if it is possible to select points (pk and ipk from the intervals Ik 
and Jk, such that |(1 + Y!k=i e^'^^'^'^'^'^^) = 0, then 

(29) max{|5-5,|} > 

Proof. Let r = maxjlS* — 5*^1}. Let us use again the "error function" 
Eix,y) = 3-1(1 + XlLi^^*''^''"'^'"'')' where Xk,yk are in 4 and 4, 
respectively. Apply for each 1 < A; < 5 the trivial estimate |e^*'^("'''~*'') — 
e2*"("'=-^'=)| < 7c{Lk + Tk) to obtain \E{x,y)\ < InJ^LiiLk + Tk) < i, 
by assumption. Also, \E{x,y)\ is a continuous function on a compact 
space, so it achieves its maximum. We claim that the maximum is 
achieved where all coordinates Xk,yk are opposite endpoints of Ik and 
Jk (i.e. if Xk is the lower endpoint of Ik then yk is the upper endpoint 
of Jk, and vice versa). Assume by contradiction that this is not so for. 
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say, Xi,yi. Then xi — yi lies in the interior of the interval /i — Ji. 
This means that for t small enough we can move Xi and/or yi to x[,y[ 
within the intervals /i, Ji so that x[ — y[ = Xi — yi + t. This yields 

(30) \E{x[,]l)\ = \S--- -e^i-i-i-y^+t) _ 1 ^ e^^'^'^-^^^l = 

k=2 

(31) = \E(x,y) + -(l-e^'''')e^'''^'='-^'^\. 

- 6 

As t varies in the neighbourhood of zero, the locus of the points E{x, y) + 
ij^g2tJ7r _ ]^-jg2i7r(a;i-yi) small arc of a circle of radius | with center 
E{x,y) — ^e^^'^^^'^~y^\ This arc goes through E{x,y) at t = 0. Com- 
bining this with the fact that \E{x,y)\ < |, it results from easy plane 
geometry that one can move along this circle in one way or the other so 
that \E{x[,y')\ becomes larger than \E{x,y)\. The same argument ap- 
plies to any of the variables Xk, yk, so we conclude that \E{x, y) \ indeed 
achieves its maximum when all Xk, yk are at some opposite endpoints of 
the intervals Ik and J^. This means that r is the maximum of \E{x,y)\. 

Therefore, the values of the function ^{^ + Ylt=i 6^*'^'-'^*="'^'=-') stay within 
a disk of radius r around S. In the first statement of the lemma we 
thus conclude that the distance of 5" from the circle of radius 
(centered at the origin) is not greater than r. In the second statement 
we conclude that the distance of S from the origin is not greater than 
r. These are equivalent to fl2^ and fl2^ . □ 

We are now ready to search for vectors u G FUB'^J'^,. For each 

fixed paii0 a~h= 1/2A^, 3/2A^, ...,{2N- l)/2N we simply take all 
possible values of ji, . . . ,75 from to A^' — 1 and check if the vector 
u given by (!^ satisfies the bound (^Hj) for all fj, (0 < j < 5). Recall 
that the coordinates of u and fj represent intervals in which the actual 
coordinates of u and fj might lie. The sum of the lengths of these 
intervals is either or + (see equations ( l26l) .(l271)). which are less 
than ^ due to our choices = 180, A^' = 19, so that Lemma [3.21 can 
indeed be applied. 

It turns out, however, that the number of vectors u satisfying the 
improved error bound (PH|) is still too high. Therefore we need to 
use the following multiscale strategy. We subdivide each interval /j^ ■* 
into two equal subintervals jj^^ ■* = [rf^^ — 1/(2A^')) ''^j^ ■*) -^],+ ^ = 
[r^^ \r^^ ^ + 1/ (2A^'))- Let rj^_ and rj^+ denote the midpoints of these 
subintervals. Clearly, each 0^ — defined in (l23l) — must fall into one 



We actually only take those a, b which approximate a, b in the triangle (0, 0) 
(1/6,0), (1/6, 1/12) i.e. in the fundamental domain. 
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of these intervals. This means that we can better approximate u by 

where e is a ± vector with the signs being chosen according to which 
subintervals (pk fall. Then is a better approximation of u and needs 
to satisfy (1281) for all fj, (0 < j < 5), with the smaller intervals 
corresponding to u^. The 2^ vectors will be called the daughters of 
u (called their mother). Clearly, if none of the daughters satisfies the 
bound fl2S|) then we can discard the mother. The point is that it often 
happens that the mother satisfies the bound fl2Hl) (corresponding to her 
own intervals), but none of her daughters do. In such a situation we 
must keep the mother at the first level of checking, but can discard 
her at the level of daughters. We then repeat this operation, obtaining 
grandchildren which have to satisfy the bound ( l28i) for all ij, (0 < 
j < 5), with even smaller intervals. Again, if none of the grandchildren 
satisfies this bound, then the grandmother is discarded. We repeat this 
operation for 7 generations, i.e. a vector u of the form survives this 
test if and only if it has a surviving descendant down to 7 generations. 
Our computer search then shows that for a fixed pair of values (a, b) 

we typically obtain 110-140 such vectors u G FUB'^J'^,. (This is quite 
satisfying considering that a, h are allowed to range over small intervals, 
and we conjecture that the precise number of unbiased vectors for any 
exact pair a, h is 48). These vectors are the candidates for the columns 

Next, we attempt to compile the basis C . If there exists a quar- 
tet {Id, F{a, b), C, D) of mutually unbiased bases, then all the columns 

of C must come from the set FUB'^^^, and, furthermore, they must 
be "almost-orthogonal" to each other. To be more precise, let u = 
■^(1, e^''^'t'\ e^^'^'t"') and V = -^(1, e^^'^'^'^, . . . , e^*'^'^^-) any two vec- 
tors from C and let 



(32) u = — (l,e ^1 ,...,e ^5 ) 
and 

(33) v= -^(l,e2*"''^i ,...,e2*"''^5^) 

be their approximation in FUB'^J'^,. Recall that we only have at hand 
the vectors u, v and keep in mind that the actual phases of u, v can lie 
anywhere in the intervals /j^ , . . . , /j^ and , . . . , /^s 5 respectively. 
The fact that (u, v) = implies the following condition on u, v: 

Definition 3.1. We will say that the vectors u and v of the form ( 132|) . 
(j33j) are A^'-orthogonal if there exist numbers (pk md ipk in the intervals 
/jf ^ and I^P, such that |(1 + J^l^ie^'^^"-^"^) = 0. 
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Similarly, we will say that u and v are N' -unbiased if there exist 
G ifj^ and e llnl\ such that + Yl=i e2-('^'=-^*))| = ^. 

These properties are clearly rotation invariant in the sense that they 
only depend on the values e ^'^^'^n ''"i , . . . ,e *'^^^^5 ''"'5 > ^ { q ^j^g 
values of rf ■* — \ • • • , t^I^ ^ —riH, ■* modulo 1. We can therefore take 

Jl ^ ' ' Jo 

mi = ■ ■ ■ = ms = and correspondingly vq = -^(1, e*'^/^', . . . , e*'^/^'), 
(where the exponents of the last 5 coordinates represent intervals, of 
course) and define the set ORT^ps^N' as the set of vectors which are 
A^'-orthogonal to vq- With this notation the rotation invariance means 
that u and v as in (l32ll , (l33ll are iV'-orthogonal if and only if the vector 

^(1, e'^'^^'n ),..., e'*"(^^5 -^-s )) is in Oi?Tep,,^,. 

Note that the set ORTeps,N' is independent of a, 6, so it can be com- 
puted once and for all at the beginning of the computer search. In or- 
der to find the set ORTf.ps^N' we first introduce the simpler set ORTj^f 
as the set of vectors u (as in flH^ ) for which there exist (pk £ Ij^. 
such that |(1 + Ylk=i^'^^'"^'') ~ 0- ^ill check each possible vec- 
tor u whether it is in ORT^' (recall that there are N'^ possibilities 
for u). In order to do so, we apply Lemma 13.21 to u and the ex- 
act vector Vq = (1, 1, . . . 1) (so that in the notations of the Lemma 
the intervals Jk are degenerate). We then use our multiscale strat- 
egy again, i.e. we test the descendants of u against vq with Lemma 
13.21 down to 7 generations. We keep only those vectors u which have 
at least one surviving descendant in each generation. Having con- 
structed the set ORTni it is now easy to obtain ORTeps,N'- Indeed, 

by definition a vector u = -^(1, e , . . . , e ^'^'^k ) can only be A^'- 
orthogonal to Vq if there exist numbers (pk in the intervals /jj, and 
ipk in [0, jij), such that \{l + YX^^e^^A'Pi.-^k)^ = 0. But then the 
numbers 0^ — ipk must fall in the intervals Ij,.-et where is either 
or 1, and hence the vector = ^(l,e ^'^^^i-^i, . . . ,e ^'^''^5-^5) is in 
ORTjyi. Therefore, to construct ORT^ps^N' we take all vectors of the 

form = :^(1, e^*'^''^i+''i , . . . , e^^'^'^^s+^s), where ek is or 1, and the vec- 



tor :^(1, e , . . . , e ^'^'"^5 ) is in ORT^i. In the specific case A^' = 19 
we found that the set ORT^ps^N' contains 322040 vectors. This means 
that the "probability" of two random vectors being A^'-orthogonal is 
322040/19^ ^ 0.13. 

Having constructed the set ORTeps^N' we search for the columns of 
the matrix C in such a way that for any two columns u and v (as 
in (l32l) . (l33l) ) we require that u, v G FUB'^^j^, and that the vector 

„. , (n') (JV'). „. , (JV') (JV') 

-ig(l, e'^'^'-n e=^*"(''« -"-"^ )) be in ORT,ps,jv 
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Let us make a last simplifying remark. It is clear that we can permute 
the columns of all appearing matrices, and hence we are free to choose 
the order of the columns of C. Therefore we assume in our search that 
the columns of C are lexicographically ordered, meaning that for any 
two columns u and v (as in ( l32l) . ( l33l) ) we have that u precedes v if 
and only if (ji, . . . , js) precedes (mi, . . . , m^) in lexicographic order. 

Our computer search has shown that for a fixed pair of values (a, h) 
there are typically 1000-5000 such A^'-orthonormal bases C . (At some 
special values of {a,b), however, there are millions. (This nicely com- 
plies with the finding of that for (a, 6) = (1/6,0) there exist 70 
bases C, whereas for generic values of (a, b) there exist only 8.) 

Finally, for any fixed pair (a, b) and any corresponding matrix C we 
attempt to compile the basis D. The columns of D must also come 
from the set FUB'^^^,, they must be A^'-orthogonal to each other, and 

they must be A^'-unbiased to the columns of C. Therefore, to find 
the candidates for the columns of D we will check any vector u e 
FUB'^\, whether it is A^'-unbiased simultaneously to all the columns 

Cfc of C. This is done by applying the trivial bound. Lemma [HTTl to the 
vector u (and its descendants for 7 generations) and the vectors for 
k = 0, . . . ,5. The reason why we use the trivial bound instead of the 
improved bound is that it speeds up calculations and very few vectors 
u survive this test anyway. Let COLf^ denote the set of surviving 
vectors, the candidates for the columns of D. If there are less than 6 
vectors in COLj^ then we conclude that D cannot exist (as it would 
need 6 columns). If there are at least 6 vectors in COLfj then we check 
whether any 6 of them can be pairwise A^'-orthogonal to each other 
(this is done by using the set ORTeps,N' again). 

Our computer search shows that there are no values of (a, b) and 
corresponding C for which all these conditions on D can be met. This 
concludes the proof of the theorem. □ 

Remark 3.3. In principle, this discretization scheme could success- 
fully be applied to settle Problem II. 1[ Of course, in the general case 
we cannot assume that B is of the form B = F{a,b). However, we 
know that B is some complex Hadamard matrix. If a complete clas- 
sification of complex Hadamard matrices of order 6 were available in 
some parametric form then a similar search could be carried out as 
above. Without such classification at hand we can still use a finite set 
of A^ representatives in each coordinate of B to approximate it with a 
quasi- Hadamard matrix B. The rest of the algorithm concerning the 
selection of quasi-unbiased vectors, and the checking of the possibly 
arising matrices C and D remains the same. Note that B has 25 free 
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entries (as the first row and column can be assumed to be 1). At first 
glance an exhaustive search for B should go through N"^^ cases, way out 
of the realm of possibilities, if ^ 100. However, one can reduce the 
number of cases with intelligent tricks, so that the search can actually 
be carried out. The problem is that while in the case of i? = F{a, b) 
and N = 180 we had only 270 candidates for B, in the general case 
we have about 10^^ such candidates already at ~ 100. And for each 
candidate B we need to run the final part of the algorithm concerning 
the selection of unbiased vectors and compiling the bases C, D. While 
it is not absolutely out of the question to carry out a computer search 
at such magnitudes, it definitely needs meticulous programming and 
probably some further mathematical ideas to reduce the number of 
cases. 

Let us also recall that the non-existence of projective planes of order 
10 was also shown by an exhaustive computer search [25] — and no 
"theoretical" proofs are known. 

4. Smooth families of mutually unbiased bases 
Throughout this section we will assume that 

(34) m = m),---m) 

is a "time-dependent" family of orthonormal bases. More precisely, 
we assume that the maps / 9 t t— > ffc(t) G (for k = 1, . . . ,d) are 
smooth (where / C M is a certain fixed open interval) and that J-'{t) 
is an orthonormal basis (ONB) for each t E I. We shall then say that 
1 1-^ {£, J-'{t)) is a smooth family of pairs of MUB if £^ = (ei, . . . e^) 
is a fixed ONB such that S and J-'{t) are mutually unbiased for all t E I . 

4.1. Smooth families of MUB and common unbiased vectors. 

Let us consider now the following question. Assume we are given a 
vector bo of unit length which is unbiased to both £ and J^(0). Can we 
"continue" bo so as to find a common unbiased vector h{t) for £ and 
JF(t) for t in a neighbourhood of 0? 

In order to answer this question we shall need some further notions 
and notations. First, if u is any vector in with ||u|| = 1, let 
be the ortho-projection on the span of u; i.e. = (x, u)u. Next, 

consider the d x d complex matrix M = [Mk^i]i<k,i<d defined by the 
formula 

(35) Mk,i := TriPkQiR) 

where Pk = Ve^, Qi = 'Pfi(o) and R = Vbo- A simple computation 
shows that 

Mfc,; = (bo,efc) (efc,f,(0))(fKO),bo). 
Note that Mk,i is unchanged if any of the vectors bo,efc,fi(0) is mul- 
tiplied by a complex number of modulus 1. Moreover, M = D1HD2 
where Di is the non-singular diagonal matrix with entries (bo,efc). 
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D2 is the non-singular diagonal matrix with entries (fi(0),bo) and 
H = [(e/j, ii{Gi))]i<k,i<d is a multiple of a Hadamard matrix. Thus M is 
of rank d. However, the real matrix N whose entries are obtained by 
taking the imaginary part of the entries of M, 

(36) Nk,i := Im(Mfe,,) = ^Tr([Pfe, Qi]R) 
cannot be of maximal rank. Indeed, for alH = 1, . . . , d, 

E = E ^^([^^' ^']^) = ^Tr([l, Qi]R) = 

since X]jt-Pfc = 1- 

Definition 4.1. Pl^e say that a unit vector b that is unbiased to both 
members of a MUB pair {S, JF) is non-degenerate if the associated real 
matrix N{h) = [lm{{h,ek){ek,ii){{i,b))]i<k,i<d has rank d - 1. 

Theorem 4.1. Let (£,J^{t)) be a smooth family of pairs of MUB. 
Let bo be a non- degenerate common normalized unbiased vector for 
the MUB pair (£,JF(0)). Then there exists an e > and a smooth 
map (—6, e) 3 t ^ h{t) such that b(0) = bo and h{t) is a common 
normalized unbiased vector to {S, J-{t)) for all t e (— e, e). 

Proof. We may modify the order of columns and rows of N :— A'"(bo) 
by rc-ordcring the vectors in S and JF. But, from the rank condition, N 
has a (d — 1) X (rf — 1) invertible submatrix so, without loss of generality, 
we may assume that the {d — 1) x {d — 1) submatrix in the upper-left 
corner of N is invertible. 

A vector v is unbiased to a given pair of MUB if and only if Av 
is so, where A G C, |A| = 1. Thus, without loss of generality, we 
may assume that (e^, bg) = rf^^/^ and in fact we may further require 
t b{t) to satisfy {ed,b{t)) = d"^/^ for all t e (-e,e). (Indeed, 
a t ^ h{t) satisfies the original requirements of our theorem, then 
t I— > h{t) :— a{t)h{t) where a{t) — d~^/'^{ed,h{t)) satisfies the just 
introduced extra condition, too.) Then the condition \{ej,h{t))\ = 
^-1/2 j = l,...,(i— 1), together with the condition of smoothness, 
are equivalent to saying that h{t) = v(a;(t)) where 

(37) v(a;) = rf-^/2|||^e-^e, + e,j 

and t i-^ x{t) e M'^"^ is a real smooth curve. In this rephrasing of 
the problem the initial condition b(0) = bo reads as x (0) = where 

= (xf\ . . . e ^'^"^ is such that {ho,ej) = d-^l'^e'^'T . Note 

also that ||v(x)|| = 1 is automatically satisfied. 
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Let us now introduce the function m : / x M'^ ^ — > M'^ ^ defined by 
the formula 

(38) Ukit, x) = I (f,(t), v(x)) I - (A; = 1, . . . , - 1) 

and note that v(x) is unbiased to J-'{t) if and only if u{t,x) = 0; that 
is, if Ukit, x) = for all k = 1, ...,d — 1. 

By assumption, u{0,x^^^) = since v(x*-°'') = bo is an unbiased vec- 
tor for T{0). Further, as |(ffc,bo)| = rf-^/^ ^ g for all k = l,..,d- 1, 
the function u is smooth in a neighbourhood of (0,a;^°'')|l A simple 
computation then shows dxi^v{x) = ■j(v(x), e^, )efc form which we de- 
duce 

dx^Uj{t,x) = |^f.(^)"^^(-3,)^| ^^ ((^XfcV(x),f;(t))(fj(t),v(x)) 

^ -Re (i{v{x),ek){ek,ij{t)){ij{t),\{x)) 



i(f.w,v(x))r 

Therefore 

dxi^Uj{t,x) = VrfRe (z(bo,efc)(efc,fj)(fj,bo) 

(t,a:)={0,x(0)) V / 

= -v^Im (^(bo,efe)(efc,fj)(fj,bo) 

Thus the Jacobian of u{0, ■) at x = x^'^^ is a nonzero multiple of the 
[d — 1) X [d — 1) submatrix in the upper- left corner of N, and the 
theorem follows by a use of the Implicit Function Theorem. □ 

Note that the Implicit Function Theorem, used in the above proof, 
actually tells us more than just existence. The invertibility of the 
Jacobian of the function m(0, ■) : — ^ ]R'^~^ at x = x^^^ guarantees 
the existence of a neighbourhood of (0, a;*-*'-*) in which the only solution 
of u{t,x) = is {t,x) = {t,x{t)). In particular, for \t\ small enough, 
in a neighbourhood of b(t), the only vectors that are unbiased for the 
MUB pair {£,T{t)) are the multiples of b(t). 

Corollary 4.2. Let h be a common non- degenerate unbiased vector for 
the MUB pair JF). Then there is a neighbourhood ofh in which all 
common unbiased vectors to [S,J-') are multiples ofh. 

Actually, our theorem also allows to prove that the number of vectors 
(counted up to multiples) that are unbiased to the family {£,J^(t)) is 
(under some non- degeneracy conditions) independent of t (for \t\ small 
enough) . 



^Actually, if s i-^ z(s) G C is smooth and z(s) ^ then — \z(s)\ 
1 

|^^Re(z'(s)z(s)). 
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Corollary 4.3. Let {S,J-'{t)) be a smooth family of pairs of MUB, and 
assume that every common normalized unbiased vector to (£^,JF(0)) is 
non- degenerate. Then there exists an e > such that for \t\ < e, 
the number of common normalized unbiased vectors to {S,J^{t)), when 
counted up to multiples, is finite and independent oft. Moreover each 



of these vectors is given by Theorem 4-1 



Proof. At t = 0, Corollary 14. 21 implies that each normalized vector that 
is unbiased to both S, and J-'{0) is isolated in the explained sense. Since 
the unit sphere of is compact, it follows that up to multiples, there 
can be only finitely many such vectors; say b^^-* , . . . , bp™'' . According 
to Theorem I4.H for each of these vectors, there is a smooth curve 
t ^ b('=)(t) such that b('=)(0) = b[,''^ and U^^t) is unbiased to {£, J^{t)). 
By the comment before Corollary 14.21 and by the fact that the just 
introduced m is a finite number, there exist some e, r > such that 
if \t\ < e then none of the vectors h^^\t), . . . , h^"^\t) are multiples of 
each other and if b is a common unbiased vector to {S, J-'{t)), then for 
every k = 1, . . . ,m, either 

||b-b(^)(t)|| > r 

or b is a multiple of h^'^^t). In particular, the number of of common 
normalized unbiased vectors to {£,J^(t)), counted up to multiples, is 
at least m (since we have the vectors h^^\t), . . . , h^^\t).) 

To prove the remaining part of our statement, we shall argue by 
contradiction. Assume there is no such e > whose existence is stated 
in our theorem. Then there should exist a real sequence t„ {n G N) 
converging to and a sequence of unit vectors b„ (n G N) such that 
for every n G N 

• b„ is a common unbiased vector to (£^,jF(t„)), 

• b„ is not a multiple of any of the vectors b*^^)(t„), . . . , h^'^\tn)- 

Since the unit sphere of is compact, it follows that there is a sub- 
sequence of hn in G N) which is convergent. In fact, without loss of 
generality we may assume that our original sequence was such. Let 
b := lim„(b„); it is then clear that ||b|| = 1 and since |(b„,efc)| = 
I (b„, fj(t„)) I = by continuity of the scalar product, absolute 

value, and the map t ^ ij{t), we have that b is a common unbiased 
vector to (£^,JF(0)). Hence by assumption there must exist a k such 
that b is a multiple of h^^^. Actually it is clear, that we even may 

(k) 

assume that b is not only a multiple of bg , but equal to it. We can 
then conclude our proof since as n — oo, we have 

r < iib^-b^'^)^*")!! ^ ||b[,')-b(')|| =0 

which is clearly a contradiction. □ 
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4.2. Unitary symmetries of mutually unbiased bases. Recall 
that if £ = (ei,...e„) and T = (fi,...f„) are two mutually unbi- 
ased ONBs, then upon multiplying each vector by a complex number 
of modulus 1 and changing the orders of the vectors in the individ- 
ual bases, they still remain two mutually unbiased ONBs. In order 
not to distinguish between such pairs, we will associate to an ONE 
S — (ei, . . . e^) a maximal abelian star algebra 

As := Span{VeJk = 1,. . . ,d} 

where Ve^ is the ortho-projection onto Ce^, {k = l,...,d). Indeed, 
As is invariant under reordering and changing phases of the vectors in 
S. Moreover, as is well known and easy to show, S = (ei, . . . e„) and 
J-" = (fi, . . . f„) are two mutually unbiased ONBs if and only if As and 
Aj^ arc quasi- orthogonal. Recall that this means that the subspaces of 
Md(C) given by 1-*- fl As and 1-*- fl Aj:- are orthogonal with respect to 
the usual Hilbert-Schmidt scalar product (^4, -B)a^^(c) = Tr {A*B) on 
Mrf(C). 

Further, to a unitary operator U on C^, we associate the authomor- 
phism au defined by the formula au{X) = UXU*. We will say that 
U implements a symmetry of {S,J-') if a{As) = As and a{Aj^) — Ajr. 
Accordingly, we shall talk about the unitary symmetry group of {£,T). 

There is a natural homomorphism from the group of symmetries of 
(£, to Sd X Sd where Sd is the group of permutations of d elements. 
Indeed, a symmetry takes a minimal projection of As and Ar into a 
minimal projection of As and Ar, respectively. Thus if a is a sym- 
metry, then there exist two permutations a — aa, /J- — /J'a & Sd such 
that 

a(^ej = i^e,(,) and a(PfJ = Pf^^,^ 
for all /c = 1, . . . ,d. Moreover, it is straightforward to show that the 
map that associates the pair {a, fj) to a, defines a group homomor- 
phism. 

Theorem 4.4. The homomorphism from the group of unitary symme- 
tries to Sn X Sn defined above is injective. In particular the group of 
unitary symmetries can have at most {d\y elements. 

Proof. For the injectivity all we need to show is that if f/ is a unitary 
operator such that UP^JJ* = Pe^ and UPf^U* = Pf^ for all A; = 1, . . . , n 
then [/ is a multiple of 1. 

However, the assumed invariance means that both the vectors of £ 
and the vectors of JF are eigenvectors for U . Thus U commutes with 
all elements of As and Aj^. As both of these are maximal abelian, it 
follows that U G As n At- However, by the quasi-orthogonality this 
intersection contains only multiples of the identity. □ 

Suppose {£,J^) is a MUB pair and that the unitary operator U im- 
plements a symmetry of {£,J^). Since U may only reorder and multiply 



22 P. JAMING, M. MATOLCSI, P. MORA, F. SZOLLOSI, AND M. WEINER 

by unit complex numbers the vectors of both 8 and JF, it is clear that if 
b is a common unbiased vector to (£, then so is Uh. Thus such uni- 
tary operators allow us to construct (possibly new) common unbiased 
vectors, once we have at least one such vector. 

Before giving a general result, let us show how this may be used to 
construct an ONB that is unbiased to both the standard basis £ and 
to the Fourier basis of C^. Let U, V be the hnear operators defined by 
the formulae 

(39) C/efc = e^^(*^-i)efc, and V^e,. = e^.i 

where the index "A; — 1" is meant by modulo d. Then f/, V are unitary, 
jjd — yd _ and VU = e^~^UV. However, by the definition of the 
Fourier basis, a simple check shows that 

(40) Vfk = e^'^^^'-'^ik, and Utk = ik+i 

for all /c = 1. . . . ,d (where the index "/c + 1" is again meant by modulo 
d). Thus U and V only change the "phases" and reorder the vectors of 
both S and J-', thus they implement unitary symmetries of {S,J-'). In 
particular, if b is a common UB vector for both the standard and the 
Fourier basis, then so is U^'V'^h for all k.l = 1, . . . ,d. 

As is well known in case of the Fourier basis, if b is a common 
normalized UB vector for then the vectors 

b, Vh, V^h, V^-^h 

form an ONB. The same stays true if one replaces V by U. 

We would like now to extend this to more general pairs of unbiased 
bases. To do so, notice first that in the above case, the natural injection 
from unitary symmetries into x Sd sends U and V to {id, a) and 
{a, id), respectively, where a & Sd is a. cyclic permutation of d, which 
is a particular example of a permutation without fixed points. 

Theorem 4.5. Let {£,J^) he a MUB pair and let h he a normalized 
vector that is unbiased to both of them. Let Uq — l,Ui, . . . ,Uk be 
unitary operators implementing symmetries cto := auo = id,ai := 
au-^, . . . ,ak '■= of {S,J-'). Assume further that for every j ^ I 
the image of aj^ o ai under the natural injection into Sn x Sn is of 
the form {a, id) or {id, a) where a & Sn is a permutation with no fixed 
points. Then 

{Uoh = b, C/ib, . . . , Ukh) 
is an orthonormal family of vectors that are unbiased to both E and T . 

Proof. Suppose that the image of aj^ o ai under the natural injection 
into Sd X Sd is of the form {a, id). Then each vector of T must be 
an eigenvector for UjUf. there exist some Xi,...,Xd G C such that 
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UlU.h = Xkh. Thus 

{Uih,Ujh) = (b,f/;f/,b) = J](b,ffc)(ffc,[/;f/,b) 

k 

= 5^(b,ffc)(f/;f/jfc,b) = 5^A,Kb,ffc)p 

k k 

For this last identity we have used the fact that b is unbiased to JF, 
thus |(b,ffc)p = 1/d, and the fact that the sum of the eigenvalues of 
a diagonalizable operator is its trace. However, by assumption U*Ui 
takes the vector into a multiple of eo-(fc); say to fik^a(^k)- Therefore 

Tr iU*Ui) = J2{U*Uiek,ek) = 5^(/ifce,(fc), e^) = 

k k 

as efc is always orthogonal to ea-[k) (since a has no fixed points). □ 

4.3. Application to the case of J^(a,b)- We shall now apply the gen- 
eral statements made so far to the case {S,J-'{a,b)) where S is the 
standard basis of C®. As it was explained, we have numerical evidence 
that up to multiple, the number of common normalized unbiased vec- 
tors is always 48. At (a, b) = (0, 0), this is a known fact. 

Theorem 4.6. There exists a neighbourhood K o/(0, 0) such that when 
counted up to multiples, {S , J-'(^a,b)) has exactly 48 common normalized 
unbiased vectors for all (a, b) G K. 

Proof. For simplicity. Theorem 14.11 and Corollary 14.31 were stated for a 
one-parameter smooth pair of MUB. However, the proofs only rely on 
the Implicit Function Theorem, so they easily extend to any number 
of parameters. 

But for {S,J-'{0,0)) i.e. for the standard and the (usual) Fourier 
basis, all 48 vectors (counted up to multiples) that are unbiased to 
them are explicitly known [7]. It is then easy (but cumbersome) to 
check that the conditions of Corollary 14.31 hold for each of them. □ 

We have seen that there is a theoretical reason (at least in a neigh- 
bourhood of the origin) behind the numerical facts that indicate that 
the number of common unbiased vectors (counted up to multiples) to 
{£, J-'{a, b)) is always 48. Unfortunately, we have so far been unable to 
find a theoretical ground for the fact that these vectors can always be 
grouped into 8 orthonormal bases. However, we may now give a partial 
result by applying what we have established about symmetries. To do 
so, first we shall need to investigate in particular the symmetries of the 
pair {£,J='(^a,b))- 

Consider the unitaries U and V defined by equation ( l39i) . For a 
generic value of the parameters a, b they do not implement symmetries. 
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However, f/^ and implement symmetries of {£, J^(a,b)) for all (a, b) G 
M?. Indeed, it is easy to check, that regardless of the value of a and b, 
we still have the relations 

U^e, = e^^2(^-i)e,, and UX = 

V^Bk = efc_3 and V% = e^^^C^-Df^^ 

where now fi, . . . , fe are the vectors of J^{a,b)- Thus by applying Theo- 
rem |42] we can draw the following conclusion. 

Corollary 4.7. Let a, b he two fixed real numbers. Suppose h is a 
common unbiased vector to the standard basis £ and to the basis J^{a,b)- 
Consider the unitaries U and V defined by equation / T^Pj) . Then all of 
the vectors in the table below 



b 













are unbiased to both bases. Moreover, each row and each column con- 
sists of pairwise orthogonal vectors. 

Unfortunately, this does not show that every common normalized 
unbiased vector can be extended to an ONB consisting of common un- 
biased vectors, only. However, in particular it shows that every com- 
mon normalized unbiased vector can be extended to an orthonormal 
triplet of unbiased vectors. 
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Appendix A. The calculations leading to Theorem 12.41 

In this section we provide the detailed calculations leading to Theo- 
rem 

Recall the form of the Fourier matrices F{0,b) from equation ([3]), 
withx= 1, ?/ = e^*'^^. We will first look for vectors u = -^(1,01,02,03,04,05). 

that obey the following further constraint^: Oi = O3O4 or equivalently 

ClCl = O3. 

We will also write 05 = r]( and 02 = fj^ with |?7| = \(\ = 1. Then 
Lemma [2.21 implies that 

Re{2c3 + v^) = 

(41) (1 + Re 03)01 + 04(1 + 03)CRer/ + (1 + c^XRer] = 

(42) (1 — Re 03)04 + ^^04(1 — 03)CIm?7 — ^^(l — 03)CIm?7 = 0. 
We will write 04 = Og and multiply (HTj) . ( l42i) by ce to obtain 

(43) (l + Reo3)ci^ + 2Re((rT^^C)Rer/ = 

(44) (1 - Re 03)c^ + 2Im ((1 - C3)ceCy)lmr] = 0. 

Notice that either fHHl) or fjl^ imply that Og is real i.e. Og = ±1. 
We will restrict our attention to Og = 1, that is Og = l,u,uj'^, thus 
04 = Og = Cg = 1, a;^ or u. 

Further, writing u = Cq(, z = —iy and using elementary computa- 
tions, we obtain 

(45) 2Rec3 + Reri^ = 

(46) (l + Reo3)(l + 2Rez/Rer7) +2Imo3lmz/Rer7 = 

(47) (1 - Reo3)(l + 2Im(z/?/) Imr]) + 2Imo3Re (z/y) Imr] = 0. 

Assume now we have a solution (03,//, 77) of the system f H5l) -( H6l) - 
( H71) . Then (1, 03O4, uffc^, 03, 04, z/t^cI) with 04 = l,a; or cu^ are solutions 
of ( |T2|) - (|T^ -( |T5|) (with x = 1). In other words, the conjugates (!) of the 
following vectors 

Wi = (l,03,Z/r/, 03,l,Z/?7), W2 = {l,C3UJ,Uf]l.j'^,C3,UJ,Ur]UJ^) 

and 

W3 = (1, 030;^ uf]uj, 03, u;^ ur]uj) 

are unbiased to both the standard and the -F(0, b) basis. It is also easy 
to see that these three vectors are orthogonal to each other. 

However, we need three more vectors. This is achieved via a miracle 
that was indicated by numerical evidence. More precisely, assume that 
(o3,r7, 1/) is a solution of the system and that {c3,fi,i>) = {—C3,iri,i)) 
is another solution of the system. We then have 6 vectors that are 
unbiased to both the standard basis and to F(0, b). Moreover, the three 

^This particular form was suggested by numerical evidence. 
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vectors Wi, W2, stemming from the first solution are orthogonal, and 
so are the ones stemming from the second solution, wl, w^, wq, namely 

W4 = (1, -C3, -zz>r/, -C3, 1, iuT]), W5 = (1, -C3U, -zz/r/u;^ -C3, u, iV7]UJ^) 

and 

Wq = {1, —Csuj'^, —ii>f]uj, —03,10'^, iu-qu). 
Finally, it is easy to check that each of 

is orthogonal to each of 

so that {wi, . . . , wq) is an orthogonal basis unbiased to both the stan- 
dard and the F(0, b) basis. 

It thus remains to exhibit two such families of solutions. To be more 
precise, we will write 77 = e** and show that, for a certain range of t, 
we may chose y = e^^^^^ in such a way that the system ( H5l) -( l46|) -(H71) 
has a solution (c3(t), e**, z^(t)), and such that there is a second solution 
{-C3{t),ie'\i>{t)). 

Now, if 1] = e**, then Rec3(t) = — 22|2l g^^^j^ g^g j^gj _ there are 

only two possibilities, csit) = ± i (l - ™^ j . For sake of 

simplicity, we will take the + sign: 

, , , , , cos2t / cos^ 
(48) C3(t) = + 



Let us first determine u = 6**^*^*^ To reduce the length and complexity 
of formulas, we will drop the dependence on t in them and simply write 

C3,/5,V5- 

But u satisfies f l46p - fH7|) which now read 

cos 2t 



- (2 — cos 2t) cos t cos (p — \/4: — cos^ 2t cos t sin (p = 1 



((2 + cos 2t) sm(3 + VA — cos^ 2t cos /?) sin t cos 

cos 2t 

+ ((2 + cos 2t) cos /? — V4 — cos^ 2t sin /?) sin t sin Lp = — 1 ^ — 

Let us write these equations in a simpler form by introducing the fol- 
lowing parameter: 



. V2 + cos 2t 
(49 j ip = arccos 
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SO that cosip = v'^+coa 2t sinip = yE^EM^ \ simple computation 
then shows that we want to solve 

(50) — sin cos t cos — cos '?/' cos t sin (y9 = — - — 

COS ip 

(51) sin(/3 + V') sintcosv? + cos(/? + 'i/') sintsin(y9 = 

Remark A.l. This system may not have solutions. For instance, it is 
easy to see that cos-^ and sin-?/' do not vanish, but the left hand side 
of (jSOll — resp. (jSTl) — vanishes when t = 7r/2 — resp. t = 0. So for t 
near or t near 7r/2, we do not expect to find a solution this way. 

The solution is now easy to obtain: 

, . — cos^ cost + cos(/5 + '?/') sin sin t 

52 cos 09 = 

^ ' ^ sin/3sin2t 

and 

. sin cos ^/^ cos t — sin(/3 + '?/') sin?/' sin t 

53 smo9 = . 

^ ' ^ sin/5sin2t 

It still has to be shown that this is a legitimate solution, that is, to 
check whether (I52l ) -( |53l) define the cosine and sine of an angle For 
this, it is sufficient to check that cos</9, sin<y9 defined by these formulas 
satisfy cos^ 99 + sin^ = 1. This easily reduces to 

cos^ cos^ t + sin^ sin^ t — 2 cos t sin t cos i\} sin il) cos /5 
(54) = sin2/?sin2 2t. 



, , . , V4 - cos2 2t 
Note that cos sm = ^ , and 

cos^ i\} cos^ t + sin^ sin^ t — sin^ 2t 

2 + cos 2t 2 2 — cos 2t 2 2 15 2 

= cos t H sin t — sin 2t = 1 — cos 2t. 

4 4 2 4 

We thus have to check that 



15 2^ V4 - cos2 2t 2 ^ 

h - cos^ 2t u + = 

2 4 4 

where u = cos (3 sin 2t. One solution of this equation is 



V4 - cos2 2t + 3 V4 - 9 cos2 2t 

cos p = 

^ 8sin2t 

V3 + sin^ 2^ + 3^9 sin^ 2t - 5 
^ ^ ~ 8sin2t 

and we omit the possible other root here. It is clear that, whenever 
^ < |sin2t| < 1 holds, we obtain a legitimate real number for 
It remains to find u = e^^ such that 

(£3(t),r7(t),z>(t)) = {-c^{t),ie'\V{t)) 
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is also a solution of 
has just been determined. 
Note that 



fH7|) . Recall, that the value of y = e 



i/3 



2Rec3 + Re77^ = -(2Rec3 + Rer^^) = 
so that is satisfied. The other two equations read 



- (2 + cos 2t) sin t cos (f + V4 — cos^ 2t sin t sin 



2 + cos 2t 



and 



( (2 — cos 2t) sin (3 — V4 — cos^ 2t cos /3) cos t cos (p+ 

2 -cos2t 



+ ((2 — cos2t) cos/? + -\/4 — cos^ 2t sin/5) cost sin;/) 

where the dependence on t in /3 and has been dropped. From this, 
we deduce that 

— sin t sin^ '?/' + cos ?/) cos t sin(/? + ?/^) 
(56) cosy? = 



(57) 



sm Lf 



sin P sin 2t 

cos i/j cos(/? + V') cos t — cos sin t sin ■j/' 



sin P sin 2t 

It is left to see that <^ is a legitimate real number, that is summing 
the squares of the two numbers defined in fl3U|) - (137|) yields 1. It is easy 
to check that this holds if and only if (1541) holds, hence there are no 
further restrictions on t. 

In summary, we have proved Theorem [2l 



Appendix B. A construction by G. Zauner that leads to 

ANOTHER ONE-PARAMETER FAMILY 



This section is inspired by G. Zauner's PhD thesis As this 

thesis is only available in German, we take this occasion to present his 
construction to a wider audience and to compare his construction to 
our family given in Theorem 12. 4[ We emphasize that the all credit for 
the results of this section goes to G. Zauner. 



Let us recall that a circulant matrix A is a matrix of the form 



( Oo 



A 



ai 
ao 



ai 



\ ai ■■■ ■ 

It is easy to check that, for each k = 0, 



ao 

0"m—l 



ai 
ao J 



2k 



{m—l)k\ 



m — 1, the vector 
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is an eigenvector of A. We may thus write A = T^ATm where Tm = 
[m~^/^ci;-'^]o<j,fc<m-i is the m x m Fourier matrix, and A is a diagonal 
matrix. 

Now, let Ao,o, ^0,1; ^1,0; ^1,1 be four m x m circulant matrices and 
write Ai 



J^^AijJ^rn where 



Aij = diag(aij(0). 



aij[Tn 



1)) : = 




is diagonal. Let T be the 2m x 2m matrix given by T 
Then 

^0,0 ^0,1 




aij{m — 1] 



T 











^1,0 A, I 









J rr 



SO that T is unitary if and only if ( ^^'"^ ) is unitary. But, this 



matrix is unitary if and only if the m matrices 5*^ 



are unitary {k = 0, ... ,m — 1). Finally, one may easily check that a 
2x2 matrix is unitary if and only if it can be written in the form 



g«/3o _j_ g«/3i 



g«/32^g«/3o _ g«/5i^ g«/32gi/33^g«/3o _|_ g*/3i^ 



For alio <k< m—1 we may thus write = S [Poi^k) , Pi{k) , P2{k) , Psi^k)) 
and define Ui = diag (e^^^(°), . . . , e^^^^"""^)) for £ = 0, 1, 2 and 3. Then 
define 

I / jr TT*-p ^ 

and 

/rq\ 771 _ fUo^m UoUsJ-'m 

^ ' V2 \U,J^n. -UiU^J'^ 

and a straightforward computation gives T = EIE2. Finally, note that 
El and E2 are Hadamard matrices so that, if T itself is a Hadamard 
matrix, then the standard matrix, the columns of Ei and the columns 
of E2 are three mutually unbiased bases in C^™. 

As an example for m = 3, Zauner [23! considers the following matrix: 



T{x) = -= 





( 1 


— e *^ 




-1 










1 


— e~"' 


ze" 


-1 




1 






1 






-1 




1 






1 










1 




— e*^ 


1 


g-JZ 








1 






1 / 
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Then T{x) is a one-parameter family of complex Hadamard matrices 



above yields a one-parameter family of MUB-triplets (/c?, Ei{x), E2{x)). 

Finally we note that Zauner's family {Id, Ei{x), E2{x)) is not equiv- 
alent to our family presented in Theorem 12.41 This can be seen in the 
following way. After dephasing the rows and columns the transition ma- 
trix T(x) = El{x)E2{x) is easily seen to be a member of the Dita-family 
D^{x) (cf. [20] for the Dita-family of complex Hadamard matrices of 
order 6). However, in our construction in Theorem 12.41 generically 
none of the appearing matrices F(0, C(t) and F(0, 6(t))*C(t) are 
members of the Dita-family. This is true, because F(0, 6(t)), C(t) are 
members of the generalized Fourier family F{a,b), while the transi- 
tion matrix F{0,b{t))*C{t) generically has a much larger Haagerup- 
invariant set than the Dita-matrices Dq{x), therefore they cannot be 
equivalent. 




Therefore, the construction 
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